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Abstract. We characterize the best L2 approximation to a multivariate function by lin- 
ear combinations of ridge functions multiplied by some fixed weight functions. In the special 
case when the weight functions are constants, we propose explicit formulas for both the best 
approximation and approximation error. 
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1. Introduction 

A function g (a- x) , where a GR n \ {0}, x £ R n , a-x is the inner product and g is a univariate 
function, is called a ridge function (in x) with the direction a. These functions and their linear 
combinations appear naturally in computerized tomography, statistics, partial differential equa- 
tions (where they are called plane waves), neural networks, and approximation theory. Ridge 
approximation in L2 was actively studied in the late 90's by K.I. Oskolkov [7], V.E. Maiorov [6], 
A. Pinkus [9], V.N. Temlyakov [10], P. Petrushev [8] and others. 

Let D be the unit disk in M?. In [5], Logan and Shepp along with other results gave a 
closed-form expression for the best L2 approximation to a function / (2:1,22) £ L2 (D) from the 
set 

ft(a\...,a r ) = j^^x) : gi :R^R, i = l,...,rj. 

Their solution requires that the directions a 1 , . . . , a r be equally-spaced and involves finite sums of 
convolutions with explicit kernels. In n dimensional case, the author [3] obtained an expression 
of simpler form for the best L2 approximation to square-integrable multivariate functions over 
some domain, provided that r = n and the directions a , ...,a r are linearly independent. 

It should be noted that problems of approximation from the set 1Z (a 1 , ...,a r ) were also con- 
sidered in the uniform norm. For example, one essential approximation method, its defects and 
advantages were discussed in [9]. Lin and Pinkus [4] characterized 1Z (a 1 , ...,a r ), i.e. they found 

r 

means of determining if a continuous function / (defined on IR n ) is of the form ^ gi (a* • x) for 

i=l 

some given a 1 , ...,a r E M ri \{0}, but unknown continuous g±, ■■■,g r - Two other characterizations 
of 1Z (a 1 , ...,a r ) may be found in Diaconis and Shahshahani [2]. Buhmann and Pinkus [1] solved 
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the inverse problem: assume that we are given a function / G TZ (a 1 , a r ) . How can we identify 
the functions gi, i = 1, ...,r? 

In this paper, we would like to consider the approximation from the more general set 

TZ (a 1 , ...,a r ; w±, ...,w r ) = j^Wj(x)^ (a 1 • x) : g { : R -> R, i = l,...,r j , 

where u; r are fixed multivariate functions. We are going to characterize the best L2 ap- 

proximation in this set (see theorem 2.4) for the case r < n. Then, in the special case when the 
weight functions wi, ...,w r are constants, we will prove two theorems on explicit formulas for the 
best approximation and the error of approximation respectively. Unfortunately, we do not yet 
know any reasonable answer to these problems in other possible cases of r. 



2. Characterization of the best approximation 

Let X be a subset of R n with a finite Lebesgue measure. Consider the approximation of a 
function / (x) = / {x\, x n ) in L2 (X) from the manifold TZ (a 1 , a r ; wi, tiv) , where r < n. 
We suppose that the functions iUj(x) and the products tUj(x) ■ </j (a* ■ x) , i = I, ...,r, belong to 
the space L2 (X) . Besides, we assume that the vectors a 1 , a r are linearly independent. We say 

r 

that a function g% = ^ w i( x )9i ( a * ' x ) m ^ ( al > ■•■■> &r ! w i> ■■■■> w r) 1S the best approximation (or 
i=i 

extremal) to / if 

ll / - 5 °llw=^ (al J a n ^ li ... >ror) ll/-^ll. 2 (x)- 

Let the system of vectors {a 1 , a r , a r+1 , a n } be a completion of the system {a 1 ,...,a r } 
to a basis in W 1 . Let J : X —> W 1 be the linear transformation given by the formulas 

yi = a l -x, i = l,...,n. (2.1) 

Since the vectors a 1 , i = 1, ...,n, are linearly independent, it is an injection. The Jacobian det J 
of this transformation is a constant different from zero. 
Let the formulas 

Xi = h l -y, i = l,...,n, 

stand for the solution of linear equations (2.1) with respect to Xj, i = 1, ...,n. 
Introduce the notation 

Y = J(X) 

and 

Yi = { yi eR: y i = a i --xi, x G X} , % = 1, n. 

For any function u £ L2 (X) , put 

n*=u* (y) = / «(b 1 -y,...,b n -y). 
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It is obvious that u* € L 2 (Y) . Besides, 

ju{y)dy = |det J\ -Ju(x)dx (2.2) 

Y X 

and 

||«*|| La(y) = |detJ| 1 / 2 .|| u || £3W . (2.3) 

Set 

Li = {w*(y)g( yi )eL 2 (Y)}, i = l,...,r. 
We need the following auxiliary lemmas. 

r 

Lemma 2.1. Zei / (x) G L2 (X). ^4 function ^ w i(^)9i ( a * " x ) ^ extremal to the function 

i=i 

r 

/ (x) if and only if ^ w*{y)g® {yi) is extremal from the space L 2 ®...®L 2 to the function f* (y). 
i=i 

Due to (2.3) the proof of this lemma is obvious. 

r 

Lemma 2.2. Let f (x) € L 2 (X). A function ^ u>j(x)g° (a* • x) is extremal to the function 

i=i 

f (x) if and only if 

J I f (x) - (a* • x) j Wj{x)h (a? • x) dx = 



i=i 



for any ridge function h (aP • x) such that Wj(x)h (a 7 • x)g L 2 (X) j = 1, ...,r. 

Lemma 2.3. T/ie following formula is valid for the error of approximation to a function 
/(x) inL 2 (X) from 1Z (a 1 , a r ; wi,...,w r ): 



E(f)=\ 11/ (x- 2 



i=i 



L 2 {X)j 



where ^ u>j(x)g? (a 1 • x) is i/ie &es£ approximation to f (x). 
i=i 

Lemmas 2.2 and 2.3 follow from the well-known facts of functional analysis that the best 
approximation of an element x in a Hilbert space H from a linear subspace Z of H must be the 
image of x via the orthogonal projection onto Z and the sum of squares of norms of orthogonal 
vectors is equal to the square of the norm of their sum. 

We say that Y is an r-set if it can be represented as Y\ x ... x Y r x Yq, where Yq is some set 
from the space M n_r . In special case, Yq may be equal to Y r+ i x ... x Y n , but it is not necessary. 
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By yW, we denote the Cartesian product of the sets Y\, Y r , Yq except for Yi, i = 1, r. That 

is, y w = Yi x ... x y^-i x y i+ i x ... x y r x y , % = 1, r. 



Theorem 2.4. Lei y 6e an r-set. A function ^2 Wj(x)<7° ( a * " x ) *s ^ e ^ es ^ approximation 

i=i 

to /(x) i/ and on/?/ i/ 



\ 



(y)dy (j) , i = i,-,r. (2.4) 



Proof. Necessity. Let a function Wi(x)gf (a* • x) be extremal to /. Then by lemma 2.1, the 



i=i 



function ^ w*(y)g < - (yi) in Ljj, © ... © L£ is extremal to /*. By lemma 2.2 and equality (2.2), 



i=i 



J r(y)w* j (y)h(y j )dy = j w*(y)h ( Vj ) £ <(y)«7? (Vi) dy (2.5) 

v v i=l 



for any product w*(y)h (yj) in L^, j = l,...,r. Applying Fubini's theorem to the integrals in 
(2.5), we obtain that 



h(yj 



J f*(y)w*(y)dyW dy 3 = j h{y 3 ) j w*(y) £ <(y) 5 ° (y^y°) 
YU) Yj yU) i=1 



dyj. 



Since h(y 3 ) is an arbitrary function such that w*Ay)h(y 3 ) G L 3 2 , 



j f*(y)w*(y)dy^= j w*(y)j2^i(y)9i(yi)dy ij \ J = l,-,r. 



Yd) 



Yd) 



Therefore, 



f*(y)-J2 w *W9i(yi 



YU) 



w*(y)dy (j) , j = l,...,r. 
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Now, since y 3 ^ Y^\ we obtain (2.4). 

Sufficiency. Note that all the equalities in the proof of the necessity can be obtained in the 
reverse order. Thus, (2.5) can be obtained from (2.4). Then by (2.2) and lemma 2.2, we finally 

r 

conclude that the function Wi{x.)gf (a* • x) is extremal to / (x). □ 

i=i 
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In the following, \Q\ will denote the Lebesgue measure of a measurable set Q. The following 
corollary is obvious. 

r 

Corollary 2.5. Let Y be an r-set. A function 9® ( a * " x ) i- n 7£ (a 1 , a r ) is the best 

i=i 

approximation to /(x) if and only if 

9 ^ y ^=\YW\j /*(y)-E^ (K)Ny 0) , j = i,...,r. 



In [3], this corollary was proven for the case r = n. 

3. Calculation of the approximation error 

In this Section, we are going to establish explicit formulas for both the best approximation 
and approximation error, provided that the weight functions are constants. In this case, since we 
vary over g^ the set 1Z (a 1 , ...,a r ; w\, ...,w r ) coincide with 1Z (a 1 , ...,a r ) . Thus, without loss of 
generality, we may assume that Wj(x) = 1 for i = 1, r. 

For brevity of the further exposition, introduce the notation 

A = J f(y)dyand f* = f*( Vi ) = j f (y)dy«, i = l,..,r. 

The following theorem is a generalization of the main result of [3] from the case r = n to the 
cases r < n. 

Theorem 3.1. Let Y be an r-set. Set the functions 

9i fa) = pi^A* " (r ~ 1) A 
and ^ 

r 

XTien £/ie function ^2 g® (a* • x) is the best approximation from 1Z (a 1 , a r ) to f (x). 
i=i 

The proof is the same as in [3]. It is sufficient to verify that the functions g® (yj) , j = 1, r, 
satisfy the conditions of corollary 2.5. This becomes obvious if note that 



r 1 1 



)l 

for j = 1, ...,r. 



/* (y)dy 



(0 



rfyO') = (r _ 1) 1 / r(y) ( /y 
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Theorem 3.2. Let Y be an r-set. Then the error of approximation to a function f(x) from 
the set TZ (a 1 , ...,a r ) can be calculated by the formula 



E( f ) = idet jp 1 / 2 (\\r \\1 2{Y) - E ^ + (r " 1} 



\y\ 



1/2 



Proof. From Eq. (2.3), lemma 2.3 and theorem 3.1, it follows that 

1/2 



s(/) = idetjr i / 2 (nrii| 2(y) -/)" 



where 



/ = 



i=l 



A_ 



L 2 (Y) 



The integral / can be written as a sum of the following four integrals: 



(3.1) 



r i - . y(i) .2 H/*llL(Y) » J 2 - |y(i)| |y(j)| j f?f] d y> 

3+i 



Y 

A 2 



It is not difficult to verify that 



/■ 



f*f*dy 



Y n x 



fc=l 



A 2 , for i,j = 1, ...,r, i / j, 



and 



k=l 



A, for i = l, r. 



(3.2) 



(3.3) 



Considering (3.2) and (3.3) in the expressions of I 2 and /3 respectively, we obtain that 



A 2 A 2 
h = r(r ~ and h = -2r(r - 1) — . 



Therefore, 



I = h + h + h + h = j2-^\\f*\\l 2{Y) -{r-i) A 



Now the last equality with (3.1) complete the proof. 



2 

in 



□ 
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Example. Consider the following set 



X = {xGM 4 :y l = y i (x)G [0;1], i = l,...,4}, 

where 

yi = xi + x 2 + x 3 - x 4 

2/2 = ari + ^2 - + £4 
y 3 = xi - x 2 + x 3 + x 4 
k 2/4 = -ari + x 2 + x 3 + x 4 



Let the function 

4 3 4 

i=l i=l j=i+l 

be given over X. Consider the approximation of this function from the set 1Z (a 1 , a 2 , a 3 ) , where 
a 1 = (1; 1; 1; — 1), a 2 = (1; 1; —1; 1), a 3 = (1; —1; 1; 1). Putting a 4 = (— 1; 1; 1; 1), we complete 
the system of vectors a 1 , a 2 , a 3 to the basis {a 1 , a 2 , a 3 , a 4 } in R . The linear transformation J 
defined by (3.4) maps the set X onto the set Y = [0; l] 4 . The inverse transformation is given by 
the formulas 

x\ = \yi + \y 2 + 32/3 - \yi 
X2 = \yi + \y2 ~ 32/3 + 32/4 
X3 = \y\ - \yi + i2/3 + 32/4 
k x 4 = -\y\ + 32/2 + 32/3 + 32/4 

It can be easily verified that /* = 2/i2/22/32/4 an d Y is a 3-set with = [0;1], i = 1,2,3. 
Besides, Yq = [0;1]. After easy calculations we obtain that A = yg; /* = |yj for i = 1,2,3; 
det J = -16; ||/*||| 2(y) = ^; Utility) = lib' * = 1 ' 2 ' 3 - Now from theorems 31 and 3 - 2 it 
follows that the function | ^=1 ( a * • x) — | is a best approximant from 72. (a 1 , a 2 , a 3 ) to / and 
/,(/) ^ V 2vl7. 
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